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Abstract 

We present a measurement of the unitarity triangle angle ^3 using a Dalitz plot analysis of the 
three-body decay of the neutral D meson from the B ± -> D^K* process. Using a 253 fb 
data sample collected by the Belle experiment, we obtain 276 signal candidates for — ► DK^ 
and 69 candidates for B^ — > D*K^, where the neutral D meson decays into Ks7r + ir~. From 
a combined maximum likelihood fit to the B ± -» DK ± and B ± -» L>*K ± modes, we obtain 
4>3 = 68° ^}jjo(stat) ± 13°(syst) ± 11° (model). The corresponding two standard deviation interval 
is 22° < <fo < 113°. 



3 



INTRODUCTION 



Determinations of the Cabbibo-Kobayashi-Maskawa (CKM) [Jj matrix elements provide 
important checks on the consistency of the Standard Model and ways to search for new 
ohysics. Various methods using CP violation in B — > DK decays have been proposed 
1, 0, 0, H[ to measure the unitarity triangle angle 3 . These methods are based on two 
tey observations: neutral D° and D° mesons can decay to a common final state, and the 
decay B + -> D^K+ can produce neutral .D mesons of both flavors via b — > cus (Fig. [T^) 
and 6 — > wcs (Fig. |IJ)) transitions, with a relative phase #+ between the two interfering 
amplitudes that is the sum, 5 + 03, of strong and weak interaction phases. For the charge 
conjugate mode, the relative phase is 6L = 5 — 03. 

(b) 1_ V "» 




B 




FIG. 1: Feynman diagrams of (a) dominant B + — > D°K + and (b) suppressed -B + — > D°K + decays 



Recently, three body final states common to D° and £)°, such as Ks^ + ^ [a], were 
suggested as promising modes for the extraction of 03. In the Wolfenstein parameterization 
of the CKM matrix elements, the amplitudes of the two diagrams that contribute to the 
decay B + -> DK + are given by M 1 ~ V* b V us ~ AA 3 (for the D°K + final state) and 
M 2 ~ V^Vcs ~ AA 3 (p + W7) (for D°K + ). The two amplitudes Mi and M 2 interfere as the 
D° and D° mesons decay into the same final state Ks7t + tt~; we denote the admixed state as 
D. Assuming no CP asymmetry in neutral D decays, the amplitude of the B + decay can 
be written as 

M+ = f(m? + , m 2 _) + re^^jim^mX) (1) 

and the corresponding amplitude for the charge conjugate B~ decay is can be written as 

M_ = f(m 2 _,m 2 + ) + re-^+^fiml, m 2 _), (2) 

where m 2 + and m 2 _ are the squared invariant masses of the Ks^ + and Ksft~ combinations, 
respectively, and /(m + ,m_) is the complex amplitude for the D° — > Ks^ + ^~ decay. The 
absolute value of the ratio between the two interfering amplitudes, r, is predicted to be 
0.1-0.2. 

Once the functional form of / is fixed by a D° — > Ksir + ir~ decay model, the D Dalitz dis- 
tributions for the B + and B~ decays can be fitted simultaneously using the above expressions 
with r, 03, and 5 as free parameters. The D° — > Ks7t + tt~ decay model can be determined 
from a large sample of flavor-tagged D° —> K s tv + tt~ decays produced in continuum e + e~ 
annihilation. 

The current measurement is an update of our analysis 0, which was based on 140 
fb _1 data sample. Both measurements are based on two decay modes, B^ — > DK ± and 
B^ — > D*K ± (D* — > Dir°). Recently, the same technique has been applied by the BABAR 
collaboration M with consistent results. 
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EVENT SELECTION 



We use a 253 fb 1 data sample, corresponding to 275 x 10 6 BB pairs, collected by the 
Belle detector. The decays B ± -> DK* and B ± -> D*K ± , D* -> Dn° are selected for the 
determination of 03; the decays B^ — > Dn^ and B^ —>■ D*n ± with D* — > Dtt° serve as 
control samples. We require the neutral D meson to decay to the Kg7r + 7r~ final state in all 
cases. We also select decays of D*^ — > Dir^ 1 produced via the e + e~ — ■> cc continuum process 
as a high-statistics sample to determine the D° — > Ks7r + 7r~ decay amplitude. The Belle 
detector is described in detail elsewhere 9| . The event selection procedures are the same as 
was described in the previous analysis [7|. 

We use AM = M K +7r - 7T + — M Ks7T + T - and Mx s7r +7r- distributions to select the — > 
Dnf events, where irf stands for "slow pion" that is a distinct signature for this decay. 
The fit to the AM distribution yields 186854 ± 856 signal events and 6126 ± 65 background 
events in the signal region (144.6 MeV/c 2 < AM < 146.4 MeV/c 2 ). The corresponding 
background fraction is 3.2%. 

The selection of B candidates is based on the CM energy difference AE = J^Ei — E heam 
and the beam-constrained B meson mass M\, c = v^bcam — (^2pi) 2 , where E heam is the CM 
beam energy, and Ei and pi are the CM energies and momenta of the B candidate decay 
products. The requirements for signal candidates are 5.272 GeV/c 2 < M\ )C < 5.288 GeV/c 2 
and \AE\ < 0.022 GeV. The AE and M bc distributions for B ± -> DK± candidates are 
shown in Fig. EJ The peak in the AE distribution at AE = 50 MeV is due to B — > 
decays, where the pion is misidentified as a kaon. The B^ 1 DK ± selection efficiency 
(11%) is determined from a Monte Carlo (MC) simulation. The number of events passing 
all selection criteria is 276. The background fraction is determined from a binned fit to 
the AE distribution, in which the signal is represented by a Gaussian distribution with 
mean AE = 0, the B^ —> Dn^ 1 component is represented by a Gaussian distribution with 
mean AE = 50 MeV and the remaining background is modeled by a linear function. The 
contributions in the signal region are found to be 209±16 signal events, 2.6±0.3 B ± -> Dt& 
events and 65±5 events in the linear background. The overall background fraction is 25±2%. 

Figure El shows the AE, M^ c and AM distributions for B ± — > D*!^ candidates. The 
selection efficiency is 6.2%. The number of events satisfying the selection criteria is 69. 
The background fraction is determined in the same way as for B ± -> events. The 

fit of the AE distribution yields 58 ± 8 signal events, 8.3 ±1.5 events corresponding to the 
linear background and 0.44 ±0.10 JB ± — > D*7r ± events in the signal region. The background 
fraction is 13 ± 2%. 



DETERMINATION OF D° -> K s iT + Tr~ DECAY MODEL 

The amplitude / of the D° —> Ks7r + 7r~ decay is represented by a coherent sum of two- 
body decay amplitudes plus one non-resonant decay amplitude, 

TV 

f(ml, m 2 _) = aje^Ajiml, m 2 ) + be ip , (3) 

3=1 

where N is the total number of resonances, ^4j(m^_, m 2 ,), a,j and ctj are the matrix element, 
amplitude and phase, respectively, of the j-th resonance, and b and (3 are the amplitude 
and phase of the non- resonant component. The total phase and amplitude are arbitrary. To 
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FIG. 2: (a) AE and (b) M^ c distributions for the — > Z) /^candidates. Dashed lines show the 
signal region. The histogram shows the data; the smooth curve in (a) is the fit result. 
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FIG. 3: (a) AE, (b) AM and (c) M bc distributions for the B ± -» D*^ candidates. Dashed lines 
show the signal region. The histogram shows the data; the smooth curve in (a) is the fit result. 



be consistent with other analyses [1,0] we have chosen the D° — > K$p mode to have unit 
amplitude and zero relative phase. The description of the matrix elements follows Ref . (TJ • 

For the D° model we use a set of 18 two-body amplitudes. These include five 
Cabibbo-allowed amplitudes: ^(892)+^-, AT*(1410) + 7r-, if *(1430)+7r-, if|(1430)+7r- and 
-fT*(1680) + 7r - ; their doubly Cabibbo-suppressed partners; and eight amplitudes with K$ 
and a vrvr resonance: K s p, K s u, K s f (980), K s f 2 (1270\ K 8 f (1370), ^(1450), K s a x 
and Ks<J 2 . The differences from our previous analysis |7( are, i) addition of if sp( 1450V, 
K *(1410) + 7r~ and its doubly Cabibbo-suppressed mode, ii) use of the Gounaris-Sakurai [12j 
amplitude description for the K$p and iCsp(1450) contributions, and iii) the mass and width 
for the / (1370) state taken from (M = 1434 MeV/c 2 , V = 173 MeV/c 2 ). 

We use an unbinned maximum likelihood technique to fit the Dalitz plot distribution to 
the model described by Eq. El We minimize the inverse logarithm of the likelihood function 
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in the form 



-21ogL 



n „ 

log p(m^_ j, m 2 J — log / p{m 2 + ,m 2 J)drn 2 + dm 

1 D 



(4) 



where i runs over all selected event candidates, and m+j, m 2 i are measured Dalitz plot 
variables. The integral in the second term accounts for the overall normalization of the 
probability density. 

The Dalitz plot density is represented by 



p(rri+,rrii) = e(rni_,rnt) I \M(rni + p 2 , nit + /f)\ 2 exp ( — ^ ) dp 2 + BOm 2 ,, m 2 

2a t mL / 



771 V 7T7T > 



(5) 

where M(m^_,m?_) = f{m 2 + ,m 2 _) is the decay amplitude described by Eq. El e{m 2 + ,m 2 _) 
is the efficiency, B{m 2 + ,m 2 _) is the background density, a m {rn^^) is the resolution of the 
squared invariant mass m? n of two pions (m 2 n = Mjj + M 2 ^ + 2M 2 — m 2 + — m 2 ). The 
free parameters of the minimization are the amplitudes aj and phases ay of the resonances 
(except for the K$p component, for which the parameters are fixed), the amplitude b and 
phase (3 of the non-resonant component and the masses and widths of the <j\ and o"2 scalars. 

The procedures for determining the background density, the efficiency, and the resolution 
of the squared invariant mass, are the same as in the previous analysis. The D° — > A's7r + 7r~ 
Dalitz plot distribution, as well as its projections with the fit results superimposed, are shown 
in Fig. HJ The fit results are given in Table |l] The parameters of the a resonances obtained 
in the fit are: M CTl = 520±15 MeV/c 2 , r CTl = 466±31 MeV/c 2 , M CT2 = 1059±6 MeV/c 2 , and 
T a2 = 59 ± 10 MeV/c 2 . The large peak in the m 2 + distribution corresponds to the dominant 
D° — > A'*(892) + 7r~ mode. The minimum in the m 2 distribution at 0.8 GeV 2 /c 4 is due to 
destructive interference with the doubly Cabibbo suppressed D° —> Ji*(892)~7r + amplitude. 
In the m 2 ^ distribution, the D° — > K$P contribution is visible around 0.5 GeV 2 /c 4 with a 
steep edge on the upper side due to interference with D° — > K$oj. The minimum around 
0.9 GeV 2 /c 4 is due to the decay D° —> K s fo(9SO) interfering destructively with other modes. 



The unbinned likelihood technique does not provide a reliable criterion for the goodness of 
fit. To check the quality of the fit, we make use of the large number of events in our sample 
and perform a binned x 2 test by dividing the Dalitz plot into square regions 0.05 x 0.05 
GeV 2 /c 4 . The test yields \ 2 = 2543 for 1106 degrees of freedom. More detailed studies are 
required in order to understand the precise dynamics of D° — > Ks^ + ^ decay. However, 
for the purpose of measuring 03, we take the fit discrepancy into account in the model 
uncertainty. 



DALITZ PLOT ANALYSIS OF B ± -> DK* DECAY 

The Dalitz plot distributions for the D — > Kstt + it~ decay are shown in Figs. 03 and |H1 for 
the B — > DK^ and B^ — > D*K ± decays, respectively. These distributions are fitted by 
minimizing the combined logarithmic likelihood function 

-21ogL = -21ogL„ -21ogL+, (6) 
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FIG. 4: (a) m?_, (b) rriL, (c) distributions and (d) Dalitz plot for the D° — > Ks7t + tt~ decay 
from the D*^ — > Dn^ process. The points with error bars show the data, the smooth curve is the 
fit result. 



where are the likelihoods of B~(B + ) data given by Eq. HJ The corresponding Dalitz 

plot densities p±(m+, w?_) are given by Eq. EJwith decay amplitudes M± described by Eq. [T] 
(B + data) and Eq. |21 {B~ data). The D° decay model / is fixed, and the free parameters of 
the fit are the amplitude ratio r and phases 03 and 5. 

We consider five sources of background (see Table Hl*|) . and determine the fraction and 
Dalitz plot shape for each component. The largest contribution comes from two kinds of 
continuum events: random combination of tracks, and correctly reconstructed neutral D 
mesons combined with random kaons. We estimate their fractions to be 21.0 ± 1.7% for 
B^ — > DK ± and 9.0 ± 2.2% for B ± — > D*K ± using an event sample in which we make 
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TABLE I: Fit results for D° — > Kstt + tt decay. Errors are statistical only. The results for the 0*1, 
o"2 masses and widths are given in the text. 



Intermediate state 


Amplitude 


Phase (°) 


Fit fraction 


K s oi 


1.57 ±0.10 


214 ±4 


9.8% 


K s p° 


1.0 (fixed) 


(fixed) 


21.6% 


K s uj 


0.0310 ±0.0010 


113.4 ± 1.9 


0.4% 


K s fo (980) 


0.394 ± 0.006 


207 ±3 


4.9% 


Ks(?2 


0.23 ±0.03 


210 ± 13 


0.6% 


K s i 2(1270) 


1.32 ±0.04 


348 ±2 


1.5% 


K 5 /o(1370) 


1.25 ±0.10 


69 ±8 


1.1% 


K 5 /j°(1450) 


0.89 ±0.07 


1 ±6 


0.4% 


K*(892)+vr- 


1.621 ±0.010 


131.7 ±0.5 


61.2% 


K*(892)-vr+ 


0.154 ±0.005 


317.7 ± 1.6 


0.55% 


A"*(1410) + vr- 


0.22 ±0.04 


120 ± 14 


0.05% 


A"*(1410)"vr + 


0.35 ±0.04 


253 ±6 


0.14% 


^(1430)+7r- 


2.15 ±0.04 


348.7 ± 1.1 


7.4% 


#5(1430)-7r + 


0.52 ±0.04 


89 ±4 


0.43% 


K 2 *(1430) + 7r~ 


1.11 ±0.03 


320.5 ± 1.8 


2.2% 


K 2 *(1430)-7r+ 


0.23 ±0.02 


263 ±7 


0.09% 


K*(1680) + vr- 


2.34 ±0.26 


110 ±5 


0.36% 


K*(1680)"vr + 


1.3 ±0.2 


87 ±11 


0.11% 


non-resonant 


3.8 ±0.3 


157 ±4 


9.7% 



TABLE II: Fractions of different background sources. 



Background source 


B ± -> DK± 


B^ -> D*K ± 


qq combinatorial 

BB events other than B ± -> D^K^/ir^ 
B± D(*) 7 r ± with K/tt misID 
Combinatorics in D° decay 
Combinatorial kaon in B^ — » D^*'K^ decay 


21.0 ± 1.7% 
2.3 ±0.2% 
0.9 ±0.1% 
0.6 ±0.1% 
<0.4% (95% CL) 


9.0 ± 2.2% 

3.1 ±0.4% 
0.7 ± 0.2% 
0.6 ±0.1% 

<0.4% (95% CL) 


Total 


25 ± 2% 


13 ± 2% 



requirements that primarily select continuum events but reject BB events. The shape of 
their Dalitz plot distribution is parameterized by a third-order polynomial in the variables 
m? + and m 2 _ for the combinatorial background component and a sum of D° and D° shapes 
for real neutral D mesons combined with random kaons. 

The background from BB events is subdivided into four categories. The D^K± and 
D^tt^ 1 combinations coming from the decay of D^*> from one B meson and K ± and 7r ± 
from the other B decay constitute the largest part of the BB background. We obtain their 
fractions of 2.3 ± 0.2% for B ± -> DK± and 3.1 ± 0.4% for B ± -> D*K± using a MC 
study. Their Dalitz plot shapes are parameterized by a second-order polynomial in the 
variables m 2 + and m 2 _ for B ± — > DK^ and by a linear function in m 2 + and m 2 _ plus D° 
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linear function in the variables m5j_ and m 2 _ plus a D° amplitude. Events in which a correctly 
reconstructed neutral D is combined with a random charged kaon are of importance. A half 
of the kaons will have the wrong sign and will be misinterpreted as decays of D mesons of 
the opposite flavor, thus introducing distortion in the most sensitive area of the Dalitz plot. 
In the MC sample, we find no events of this kind, which allows us to set an upper limit of 
0.4% (at 95% CL) on the fraction for this contribution. 

To test the consistency of the fitting procedure, the same fitting procedure is applied 
to the Br* 1 — > control samples as to the B ± — > D^K ± signal. In the case of 

B± dWtt*, a small amplitude ratio is expected (r ~ \V u bV* d \/\V cb V* d \ ~ 0.01 - 0.02). 
Here, we consider B + and B~ data separately, to check for the absence of CP violation. 
The free parameters of the Dalitz plot fit are r± and 9±, where 9± = 5 ± 03. The fit results 
for B ± -> Dir* sample (3425 events) are r + = 0.039 ± 0.021, 9 + = 240° ± 28° for B + data 
and r_ = 0.047 ± 0.018, 6L = 193° ± 24° for B~ data. It should be noted that, since the 
value of r is positive definite, the error of this parameter does not serve as a good measure 
of the r ~ hypothesis. To demonstrate the deviation of the amplitude ratio r from zero, 
the real and imaginary parts of the complex amplitude ratio re 10 are more suitable. Figure [7| 
(a) shows the complex amplitude ratio constraints for the B + and B~ data separately. It 
can be seen that the amplitude ratios differ from zero by roughly two standard deviations 
for both cases but are not inconsistent with the expected value of r ~ 0.01 — 0.02. 

The other control sample, B — > D*^ with D* decaying to Dtt , also does not show 
any significant deviation from r ~ 0. The results of the fit to the B^ — > D*^ sample (642 
events) are r + = 0.015 ± 0.042, 9 + = 169° ± 186°, r_ = 0.086 ± 0.049, 9_ = 280° ± 30° and 
are shown in Fig. [7| (b). 

RESULTS 

Fig. IS] shows the constraints on the complex amplitude ratio re %e for the B ± — > DK ± 
and B ± — > D*K ± samples. It can be seen that in both signal samples a significant non-zero 
value of r is observed. A difference between the phases 9 + and 6L is also apparent in both 
the B^ — > DK^ and B^ D*K ± samples, which indicates a deviation of 03 from zero. 

A combined unbinned maximum likelihood fit to the B + and B~ samples with r, 03 and S 
as free parametejs yields the following values: r = 0.25 ±0.07, 03 = 64° ±15°, 5 = 157° ±16° 
for the 5 ± -> sample and r = 0.25 ± 0.12, 3 = 75° ± 25°, 5 = 321° ± 25° for the 

B^ — > D*K ± sample. The errors quoted here are obtained from the likelihood fit. These 
errors are a good representation of the statistical uncertainties for a Gaussian likelihood 
distribution, however in our case the distributions are highly non-Gaussian. In addition, the 
errors for the strong and weak phases depend on the values of the amplitude ratio r {e.g. 
for r = there is no sensitivity to the phases). A more reliable estimate of the statistical 
uncertainties is obtained using a large number of MC pseudo-experiments as discussed below. 

The model used for the D° — > A's7r + 7r~ decay is one of the main sources of systematic 
error for our analysis. Since the Dalitz density is proportional to the absolute value squared 
of the decay amplitude, the phase (p{m\,m 2 _) of the complex amplitude is not directly 
measured. The phase variations across the Dalitz plot distribution are therefore the result 
of model assumptions. 

A MC simulation is used for estimating the effects of the model uncertainties. Event 
samples are generated according to the Dalitz plot distribution described by the amplitude 
given by Eq. Qwith the resonance parameters extracted from our fit of continuum D° data. 
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FIG. 7: Constraint plots of the complex amplitude ratio re ld for (a) B — > Dir^ and (b) — > 
L)*7r ± decays. Contours indicate integer multiples of the standard deviation. Dotted contours are 
from B + data, dashed contours are from B~ data. 
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FIG. 8: Constraint plots of the complex amplitude ratio re 1 ® for (a) B^ — > DK^ and (b) B^ — » 
D*K^ decays. Contours indicate integer multiples of the standard deviation. Dotted contours are 
from B + data, dashed contours are from B~~ data. 

We then fit the distribution using different models for f(m + ,m_) (see Table ITTTj) . We scan 
the phases 03 and 6 in their physical regions and take the maximum deviations of the fit 
parameters ((Ar) max , (A0 3 ) max , and (A<5) max ) as model uncertainty estimates. The values 
for (Ar) max , (A0 3 ) max and (A<5) max quoted in Table ITTT1 are obtained with the value r = 0.13. 
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TABLE III: Estimation of model uncertainty. 



Fit model 


(Ar) max 


(A<fe)max (°) 


(A5) max (°) 


F r = F D = l 


0.01 


3.1 


3.3 


T(q 2 ) = Const 


0.02 


4.7 


9.0 


Narrow resonances plus non-resonant term 


0.03 


9.9 


18.2 


Total 


0.04 


11 


21 



For larger r values, the model uncertainty tends to be smaller, so our estimate of the model 
uncertainty is conservative. 

All the fit models are based on Breit-Wigner parameterizations of resonances. Since a 
Breit-Wigner amplitude can only describe narrow resonances well, the usual technique to 
deal with broad states is to introduce Blatt-Weisskopf form factors for the D° meson (Fd) 
and intermediate resonance (F r ) and a g 2 -dependence of the resonance width T. These 
quantities have substantial theoretical uncertainties. We therefore use a fit without Blatt- 
Weisskopf form factors and with a constant resonance width to estimate such an error. We 
also use a model containing only narrow resonances (if* (892), p, doubly Cabibbo-suppressed 
if* (892) and /°(980)) with the wide ones approximated by the flat non-resonant term. The 
study of the model errors is summarized in Table IIII1 Our estimate of the systematic 
uncertainty on 3 is 11°. 

There are other potential sources of systematic error such as uncertainties in the back- 
ground Dalitz density, efficiency variations over the phase space, resolution, and possible 
fit biases. These are listed in Table ITVl for the B ± — > DK ± and B^ — > D*K ± modes sepa- 
rately. The effect of background Dalitz density is estimated by extracting the background 
shape parameters from the M& sidebands and by using a flat background distribution. The 
maximum deviation of the fit parameters from the "standard" background parameterization 
is assigned as the corresponding systematic error. The effect of the uncertainty in the back- 
ground fraction is studied by varying the background fraction by one standard deviation. 
The efficiency shape and resolution are extracted from the MC simulations. To estimate 
their contributions to the systematic error, we repeat the fit using a flat efficiency and a fit 
model that does not take the resolution into account, respectively. The biases due to the 
efficiency shape differ for — > DK ± and B — > D*K ± samples, but since we expect the 
values of the efficiency systematics to be close for the two modes, we assign the maximum 
value of the bias as the corresponding systematic error. 

We use a frequentist technique to evaluate the statistical significance of the measurements. 
To obtain the probability density function (PDF) of the fitted parameters as a function of 
the true parameters, which is needed for this method, we employ a "toy" MC technique that 
uses a simplified MC simulation of the experiment which incorporates the same efficiencies, 
resolution and backgrounds as used in the data fit. This MC is used to generate several 
hundred experiments for a given set of r, 9 + and 6L values. For each simulated experiment, 
Dalitz plot distributions are generated with equal numbers of events as in the data, 137 
and 139 events for B~ and B + decays, correspondingly, for B ± — > DK^ mode and 34 and 
35 events for B~ and B + for B^ — > D*K ± mode. The simulated Dalitz plot distributions 
are subjected to the same fitting procedure that is applied to the data. This is repeated 
for different values of r, producing distributions of the fitted parameters that are used to 
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TABLE IV: Contributions to the experimental systematic error. 



Source 


B ± -> DK* 


B 


± -> D*K^ 


Ar 


A0 3 (°) 


as n 


Ar 


A0 3 (°) 


AS O 


Background shape 


0.027 


5.7 


4.1 


0.014 


3.1 


5.3 


Background fraction 


0.006 


0.2 


1.0 


0.005 


0.7 


1.4 


Efficiency shape 


0.012 


4.9 


2.4 


0.002 


3.5 


1.0 


m 2 ^ resolution 


0.002 


0.3 


0.3 


0.002 


1.7 


1.4 


Control sample bias 


0.004 


10.2 


10.2 


0.004 


9.9 


9.9 


Total 


0.03 


13 


11 


0.02 


11 


11 



produce a functional form of the PDFs of the reconstructed values for any set of input 
parameters. 

We parameterize the PDF of a set of fit parameters (r, 03,5), assuming the errors of 
parameters R,e(r±e l6± ) and lm(r±e ld± ) are uncorrelated and have Gaussian distributions 
with equal RMS which we denote as a. The PDF of the parameters (r±, 0±) for the true 
parameters (f±, 9±) is thus written as 



d 2 P(r±,9±\r±,9±) = 7j^ ex P 



(r ± cos 9± — r cos 9± ) 2 + (r± sin #± — r sin #± ) 2 1 
— r±dr±d9±. 

(7) 

To obtain the PDF for the parameters (r, 3 , S) we fix r = r + = r_ and substitute the total 
phases with <5 + 3 and 5 — 03 : 

cf 3 P --(i 2 P - - d 2 P 

(r, 03, <S|f, 03, 5) = - — — (r, 5 + 3 |f, 5 + 3 )- — — — (r, 5 - 3 |f, 5 - 3 ). (8) 



drdfodS dr + d9 + dr_d9^ 

There is only one free parameter a which is obtained from the unbinned maximum likelihood 
fit of the MC distribution to Eq. |HJ The value of a is equal to 0.10 for the B ± -> 
decay and a = 0.18 for the B ± — » D*K ± decay. 

Once the PDF is obtained, we can calculate the confidence level a for each set of true 
parameters (f, 3 , 8) for given measurements, (r, 3 , 5) = (0.25, 64°, 157°) for the DK^ 1 mode 
and (0.25, 75°, 321°) for the D*K ± mode. The confidence regions for the pairs of parameters 
(03, S) and (03, r) are shown in Fig.H( J B ± -> DI^ mode) and Fig. EH (P^ -> D*K^ mode). 
They are the projections of the corresponding confidence regions in the three-dimensional 
parameter space. We show the 20%, 74% and 97% confidence level regions, which correspond 
to one, two, and three standard deviations for a three-dimensional Gaussian distribution. 

For the final results, we use the central values that are obtained by maximizing the 
PDF and the statistical errors corresponding to the 20% confidence region (one standard 
deviation). Of the two possible solutions (03, 5 and 03 + 180°, 5 + 180°) we choose the one 
with < 03 < 180°. The final results are 

r = 0.21 ± 0.08 ± 0.03 ± 0.04, 3 = 64° ± 19° ± 13° ±11°, 5= 157° ± 19° ± 11° ± 21° (9) 

for the JB ± — > DK ± mode and 

r = 0.12+°;^ ± 0.02 ± 0.04, 3 = 75° ± 57° ± 11° ± 11°, 5 = 321° ± 57° ± 11° ± 21° (10) 
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for the — > D*K mode. The first, second, and third errors are statistical, systematic, 
and model dependent errors. 

While the 03 and 5 values that are determined from the toy MC are consistent with those 
that are determined in the unbinned maximum likelihood fits for both DK* and D*K±, 
the corresponding r values are significantly different. This is due to a bias in the unbinned 
maximum likelihood. Since r is a positive-definite quantity, the fit tends to return a larger 
value for r than its true value, particularly when r is small. 

In the frequentist approach, the significance of CP violation is evaluated by finding the 
confidence level for the most probable CP conserving point, i.e. the point with r = or 
03 = 0, for which the confidence level a(f, 03, 5) is minimal. This procedure gives a = 94% 
for the B ± -> DK± sample and a = 38% for B ± -»■ D'K*. 

The two events samples, 5 ± — > DK^ and B ± — > D*K ± , are combined in order to 
obtain a more accurate measurement of 03. The technique we use to obtain the combined 
measurement is also based on a frequentist approach. Here we have five true parameters (03, 
Ti, Si and #2, where the indices 1 and 2 correspond to B^ 1 — > DK^ and B — > D*K ± 
modes, respectively) and six reconstructed parameters (r, 03 and 5 for each of the two 
modes). The PDF for the reconstructed parameters is written as 

-J-( X ^) = —FlT^Jr( r ^ (03)l, $l\fl, 03, Si) d ^ B ,~2 D *l K (r 2 , (03)2,52^2, 03, £2), (H) 

ax dracpzdo drdcp^do 

where x = {dr X) d(4> 3 )i, d6i, dr 2 , c?(0 3 )2, dS 2 ) is a vector of the reconstructed parameters, and 
fi = (03, fx, f 2 , S\, 5 2 ) is a vector of the true parameters. Using this PDF and Feldman- 
Cousins likelihood ratio ordering 0], we can calculate the confidence level a(fi). This 
approach gives 03 = 68° as the central value. The one standard deviation interval for 03 
(which corresponds to the 3.7% confidence level for the case of a five-dimensional Gaussian 
distribution) is 03 = 68° l 15 o. 

Since the B^ — > DK^ contribution dominates in the combined measurement, we use its 
value of the systematic uncertainty, which is 13°, as an estimate of the systematic uncertainty 
in the combined 03 measurement. The 03 result from the combined analysis is 

03 = 68° t\£o ± 13° ± 11°, (12) 

where the first error is statistical, the second is experimental systematics, and the third 
is model uncertainty. The two standard deviation interval including the systematic and 
model uncertainties is 22° < 03 < 113°. The statistical significance of CP violation for the 
combined measurement is 98%. 



CONCLUSION 

We report results of a measurement of the unitarity triangle angle 03 that uses a method 
based on a Dalitz plot analysis of the three-body D decay in the process 
The measurement of 03 using this technique was performed based on 253 fb^ 1 data sample 
collected by the Belle detector. From the combination of B^ — > DK^ and B^ — >• D*K ± 
modes, we obtain the value of 3 = 68° ±{£ ± 13° ± 11° (solution with < 3 < 180°). 
The first error is statistical, the second is experimental systematics and the third is model 
uncertainty. The two standard deviation interval (including model and systematic uncer- 
tainties) is 22° < 03 < 113°. The statistical significance of CP violation for the combined 
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FIG. 9: Confidence regions for the pairs of parameters (a) (r, 4>s) and (b) (^3, 5) for the — > DK^ 
sample. 
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FIG. 10: Confidence regions for the pairs of parameters (a) (r, ^3) and (b) (4>3,S) for the 
D*!^^ sample. 



measurement is 98%. The method allows us to obtain a value of the amplitude ratio r, 
which can be used in other 03 measurements. We obtain r = 0.21 ± 0.08 ± 0.03 ± 0.04 for 
the B ± -> DK± mode and r = 0.12j£u ± 0.02 ± 0.04 for the B ± -> D*K± mode. 
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